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Abstract 

This paper considers random Hill's equations in the limit where the periodic 
forcing function becomes a Dirac delta function. For this class of equations, the 
forcing strength the oscillation frequency and the period (Ar)^ are allowed 
to vary from cycle to cycle. Such equations arise in astrophysical orbital problems in 
extended mass distributions, in the reheating problem for inflationary cosmologies, 
and in periodic Schrodinger equations. The growth rates for solutions to the periodic 
differential equation can be described by a matrix transformation, where the matrix 
elements vary from cycle to cycle. Working in the delta function limit, this paper 
addresses several coupled issues: We find the growth rates for the 2x2 matrices that 
describe the solutions. This analysis is carried out in the limiting regimes of both large 
qk S> 1 and small qt <C 1 forcing strength parameters. For the latter case, we present 
an alternate treatment of the dynamics in terms of a Fokker-Planck equation, which 
allows for a comparison of the two approaches. Finally, we elucidate the relationship 
between the fundamental parameters (Xk,qk) appearing in the stochastic differential 
equation and the matrix elements that specify the corresponding discrete map. This 
work provides analytic — and accurate — expressions for the growth rates of these 
stochastic differential equations in both the q^ 3> 1 and the q^ <C 1 limits. 
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I. INTRODUCTION 



This paper considers random Hill's equations in the delta function limit. A random 
Hill's equation can be written in the form 

^| + [A fc + g k Q(t)]y = , (1) 

where the barrier shape function Q(t) is periodic, so that Q(t + At) = Q(t), where At is 
the period (we generally take At = ir). In the delta function limit, the periodic functions 
Q(t) become Dirac delta functions, 

Q(t) = 5([t] - tt/2) , (2) 

where the square brackets indicate that time is measured mod-7r. The parameter q k denotes 
the forcing strength, which is a random variable that takes on a new value every cycle 
(where the index k determines the cycle). The parameter Afc, which determines the natural 
oscillation frequency of the system, also varies from cycle to cycle. In this work the period 
At is considered fixed; one can show that cycle to cycle variations in At can be scaled out 
of the problem and included in the distributions of the (A*., qk) [2]. 

The original form of Hill's equation (OQ) holds the values of the parameters constant 
[13], and such equations arise often in physics [19]. A straightforward generalization of 
this classic problem is to consider parameters that vary from cycle to cycle (according to 
a well-defined distribution). Further, the limit of delta function barriers (equation [2]) 
arises in many applications (see below) and thus provides a natural starting point for this 
analysis. 

One specific motivation for considering random Hill's equations arises from orbit prob- 
lems in astrophysical settings, including dark matter halos, galactic bulges, tidal streams, 
and young embedded star clusters. These astrophysical systems generally have non- 
spherical, extended mass distributions, with corresponding potentials that are asymmetric. 
With this loss of symmetry, angular momentum is not conserved, orbits are not confined 
to particular planes, and orbital instabilities often arise. For example, if an orbit is ini- 
tially confined to the principal plane of a dark matter halo (or any triaxial, extended mass 
distribution), the motion is unstable to perturbations out of the orbital plane (see Ref. [4] 
and Appendix A). The development of this instability [2] is described by a random Hill's 
equation (as given by equation [Tj), with sharply-peaked forcing barriers that can be de- 
scribed by delta functions (as given by equation [2]). This orbit instability arises in many 
other astrophysical systems, including embedded young star clusters, galactic bulges, and 
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tidal streams [4,6]. This instability produces a number of astrophysical effects, including 
changing the velocity distributions from radial to more isotropic, making highly flattened 
systems more rounded, and helping to disperse tidal streams. 

In another application, the reheating epoch at the end of the inflationary phase in the 
early universe [12] is described by a "parametric resonance instability" [17]. During infla- 
tion, the potential of the inflaton field $ dominates the energy density, which is primarily 
in the form of vacuum energy. After sufficient inflation has taken place, this energy must be 
converted into matter and radiation so that the universe can evolve into its present state. 
As a result, the inflaton field $ must couple to matter fields \i an d the subsequent conver- 
sion of energy is governed by a Hill's equation, often the Mathieu equation [15]. Additional 
fluctuations [15,20] that are present during this process — due to thermal and quantum 
effects — convert the equation of motion for reheating into a random Hill's equation (see 
Appendix B). During inflationary reheating, the fluctuations have small amplitudes and 
the forcing terms can be modeled as delta functions, so that the resulting problem is de- 
scribed by equations (CQ) and (j2j). The instability itself acts to rapidly convert the vacuum 
energy of the universe into matter and radiation. In the classical problem, however, the 
parameter space of Hill's equation retains bands of stability that can inhibit this conver- 
sion. Random fluctuations tend to erase these bands of stability, as shown herein, and 
thereby increase the efficacy of the reheating process. 

For completeness, we note that in quantum systems with periodic lattices and a source 
of noise, the corresponding Schrodinger equation takes the form of a random Hill's equation 
[5,7]. This topic is relatively well developed [20], but the results presented in this paper 
can also be useful in this context. 

Periodic differential equations in this class can be described by a discrete mapping of 
the coefficients of the principal solutions from one cycle to the next. The transformation 
matrix takes the form 

_ \h k (h\ - l)/g k ] 
Mk -[g k h k J' ^ 

where the subscript denotes the cycle. The matrix elements for the kth cycle are given by 



h k = yi(n) and g k = y^n) , (4) 

where y\ and yi are the principal solutions for that cycle. For Hill's equation with delta 
function barriers (121) . the principal solutions have been found previously [4] and the matrix 
elements take the form 

h k = cosy? fc - q t- sin and g k = -Jx^sin(p k - q k cos 2 ((p k /2) , (5) 
2v \ k 
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where we have defined ifk = \/Afc7T. The index k indicates that the quantities (Afc, qk), and 
hence the solutions (hk,gk), vary from cycle to cycle. Throughout this work, the random 
variables are taken to be independent and identically distributed (iid). 

Note that one might expect the matrix in equation ([3]) to have four independent ele- 
ments (instead of only two). In this paper, however, we specialize to the case where the 
periodic functions Q(t) are delta functions, which are symmetric about the midpoint of the 
period. This property implies that 2/1(71") = y<i (it), which eliminates one independent matrix 
element [1,19]. In addition, since the Wronskian of the original differential equation (pQ) 
is unity, the determinant of the matrix map must be unity, and this constraint eliminates 
another independent element. 

The growth rates for Hill's equation (JTJ are determined by the growth rates for matrix 
multiplication of the matrices M.^ given by equation Q. Here we denote the product of 
N such matrices as M> N \ and the growth rate 7 is defined by 

7= Km ^\og\\M^\\. (6) 
Previous work [10,11,18] shows that this result is independent of the choice of the norm 

This paper is organized as follows. To determine the growth rates for the differential 
equation (JTJ), we consider the multiplication of infinite strings of random matrices of the 
form ([3]). Working in the delta function limit, this paper explores the regime where the 
forcing strengths are large qk 3> 1 (Section II), and the opposite regime of small forcing 
parameters qk <C 1 (Section III). For the small qk limit, we develop an alternate treatment 
of the dynamics using the Fokker-Planck equation (Section IV). Next, we consider the 
relationship between the matrix elements {hk,gk) that appear in the discrete map ([3]) and 
the random variables (Xk,qk) that appear in the original differential equation (pQ). For the 
limiting case of delta functions barriers, we find this transformation explicitly (see equation 
[5]), and constrain the distributions of the matrix elements for given distributions of the 
input parameters in Section V. The paper concludes (in Section VI) with a summary of 
the results and a brief discussion of future applications. In addition to the appendices 
that outline the physical motivation for random Hill's equations, we also present a simple 
iterative map (Appendix C) that reproduces our basic results for the growth rates. 

II. THE LIMIT OF LARGE FORCING STRENGTH PARAMETER 

This section considers the case of large forcing strengths qk for the problem with delta 
function barriers. This limit applies to Hill's equations that govern the orbit instabilities 
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in dark matter halos, embedded young star clusters, and other extended mass distributions 
(Appendix A). For the triaxial orbits [4] that originally motivated this study, for example, 
the forcing strengths q k ~ 1000 when the period Ar^ and oscillation parameter are of 
order unity. 

In this limit, it is useful to factor the matrix M. k that connects solutions from cycle to 
cycle by writing it in the form 



M. k = h k B k where Bk 



l/x k 1 



(7) 



where x k = h k /g k and where 4> k = 1 — \/h\. The ansatz of equation ([7j) separates the 
growth rate for this problem into two parts: 7 = 7^ + 7b- The first part 7^ of the growth 
rate is given by 

1 N 

lh = lim — £>g|/i fc | , (8) 
N ^°° N k =i 

while the remaining part 7b is determined by matrix multiplication of the matrices B k . In 
the general problem [2,3], much of the work thus involves finding the growth rates 75. In 
the limit of large q k , however, much simpler — but approximate — forms for the growth 
rates can be found, as shown below. 

Theorem 2.1: Consider a random Hill's equation with a delta function barrier. In the 
limit of large q k 3> 1 and constant A, the growth rate has the form 

7 =<log|2/i fc |) +0(l/q k ), (9) 

where the h k are given by equation (jlj) and where the angular brackets denote the expec- 
tation value. This form is valid provided that A 7^ n 2 , where n is an integer. The width of 
the zone for which equation is not valid has order 5X = 0(l/q k ). 
Proof: We separate the problem into two pieces according to equation The quantities 
that appear in the matrix elements in the B k are the ratios x k = h k /g k and the correction 
factors (j) k = 1 — l/h k . For the case of delta function barriers, considered here, the quantities 
Xk and (p k can be written in the form 

q k (n/<p) siny? -2 cosy) / %p \ 2 

x k = — 7T- r— T7 1 v . and fa = 1 - : , 10 

q k (l + cos^jj + 2{(p/7r) sm(p \^qk sm(p — zcpcoscp J 

where we have suppressed the subscripts on the angles (p = y/Xn. In the limit of large 
forcing strength q k , the ratios x k become independent of the values of the q k . In particular, 
x k and (p k take the asymptotic forms 

7rsin(/9 

lim x k = — -. r and lim m k = 1 . (11) 

Qk^oo ip[\ -f COS(^) Qk-^oo 
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For constant A, the angles ip = a/A~7t are also constant, and the x k are all the same in the 
limit q k — > oo. In this same limit, the matrices B k become constant from cycle to cycle 
and are denoted here as B . These matrices have the simple multiplication property 

B 2 = 2B , (12) 

so that the growth rate 7# = In 2. The remaining part of the growth rate 7^ is given by 
equation (jSJ), where the h k are given by equation Combining these two results yields 
the expression for the growth rate given in equation 0. The correction term is considered 
below. 

The above derivation of the growth rate 7^ is valid as long as the x k in equation ffTTl) 
remain finite and nonzero. This requirement leads to the condition that A ^ n 2 , where n 
is an integer. This result also follows from a previous theorem of Ishii [14,20]. We thus 
obtain the stated restriction on the range of validity of the growth rate. To constrain the 
width of the angular zone for which this result is not valid, we write ip = mr + Sip. Finding 
the condition for which \h k \ < 1, which is the condition for one cycle to be stable, we find 
that 

Sip < = — . (13) 

The width of the angular zone thus depends on q k , which varies from cycle to cycle, but 
Sip = 0(1/ qk)- The width of the stability zone for A is given by 

S\ = — Sip = . (14) 

vr irq k 

The width of the zone (in the parameter A) for which the growth rate of equation (Q fails 
is thus of order l/qk, as claimed. 

In deriving the leading order term in equation (jH]), we have considered the variables x k 
to be constant in the limit of interest. In the more general case, the Xk vary from cycle 
to cycle (note that varying X k , not considered here, would also contribute to variations in 
x k ). Including these variations leads to a correction A7 to the growth rate. As shown in 
Theorem 2 of Ref . [2] , this correction can be written in the form 

A 7 = (\og\l + x kl /x k2 \) -log 2, (15) 

where the x k i and x k 2 represent two independent samples of the variable x k (see equation 
[TO]). Note that A7 — > in the limit where the x k are constant. In the limit where the 
q k are large, but not infinite, variations in the x k are small, and equation f|T5l) can be 
expanded and written in the approximate form 

^kl <Pk2 \ , ^, _2 




q k i sm ip k i q k2 sm ip k2 



+ 0(q k 2 ) - (16) 
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The first term is 0(1/ qk)- If the distributions of q k and (p k are symmetric, the first term 
can vanish, and the correction A7 to the growth rate becomes second order in l/qfc. 

The result given in equation (TT5I) includes the variations of the x k but does not take 
into account possible deviations of the <fi k from unity. In order to determine how large 
these corrections can be, we consider the case where the correction factors <fi k are close to 
- but not exactly - unity If 7^ is the true growth rate for multiplication of the matrices 
Bk, and 70 is the growth rate obtained in the limit where — > 1, then we denote the 
difference as £7 = 70 — 75. For delta function barriers, we can use Theorem 2.3 of Ref. [3] 
to write this correction term in the form 

1 N r 2 4-tn 2 1 

5 7 = hm i- V % 1 - , (17) 

N ^[ (xk + Xk+i ) [Xk + Xk-i ) sin' ip k i\ 2 q{ 

where (p k — V^fc 71 "- This expression is correct to leading order in l/q k . As a result, we can 
use the asymptotic expressions for the Xk to evaluate £7 when the q k are large (but not 
infinite). In this case, the x k are independent of the qk, and the correction £7 has the order 



5 7 = O (\ k /q 



where <p k = yAfc 71 "- As the q k become large, this correction decreases, and we recover the 
growth rates given by equations (El) and ( fi~5l) . If the angle <p is held fixed, then the x k in 
equation ffTTj) are identical in the large q k limit, and this correction £7 to the growth rate 
reduces to the form 

^(^y 2 (L)^/L\. (19) 

Y 71 sm ip J \q 2 k / sin> \q 2 k / 

Note that these results are only valid when sirup 7^ 0, which in turn requires A^n 2 (where 
n is an integer). • 

The above considerations provide corrections for the growth rate of Theorem 2.1 for 
cases where the q k are large, but corrections of order l/q k are still relevant. In the opposite 
limit where q k — > 00, the growth rate can be simplified further: 

Corollary 2.1: For a random Hill's equation with delta function barriers, in the limit 
q k — > 00, the growth rate approaches the form 



7 = ( log 



qk . rr 
—j= sm v A7T 



(20) 



This form is valid for A/n 2 , where n is an integer. 

Proof: This limit represents a stronger version of the conditions for which equation ([9j) is 
valid. Starting with equation (jlj), the matrix elements h k approach the following form in 
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the limit g& — > oo: 

hk — > — -4= sin 09. (21) 

Using this expression in equation ([9]), we obtain the claimed form for the growth rate given 
in equation ( J201) . • 

Notice that the difference between the asymptotic form of the growth rate from equation 
f T2"0"j) and the approximation of equation is first order in l/g&. For comparison, the 
corrections due to variations in the Xk are first order for asymmetric variations and second 
order for the symmetric case (see equation [16] ). The corrections due to the departure of 
the 0fc from unity are second order (equations [UJ and [IB]). 

Figure [H illustrates the validity of the approximations derived in Theorem 2.1 and its 
corollary. Here we take Hill's equation in the delta function limit, with fixed oscillation 
parameter A. We then let the forcing strength vary according to = (1 + £,k)%, where 
is a uniformly distributed random variable and the constant go determines the amplitude 
of the fluctuations. The growth rates are calculated using three successive approximations: 
(a) the full matrix multiplication scheme from equation ([3]) with growth rate given by 
equation ([6]), (b) the approximation of equation (Q which is valid in the limit of large 
forcing strength g&, and (c) the more extreme approximation of equation ( l20i) . which is 
valid in the limit qk — > oo. The upper curve in the figure shows the difference between 
the full growth rate of equation ([6]) and that of equation (TSUI) ; the lower curve shows the 
difference between the full growth rate and that of equation (J9]). Both approximations work 
well for large qk, measured here using qo, where "large" means go greater than ~ 100. Both 
curves approach power-law forms, with well-defined slopes, showing that the approximation 
of Theorem 2.1 is valid to second order in 1/qo, whereas the more extreme approximation 
of Corollary 2.1 is only accurate to first order in 1/go- 

III. THE LIMIT OF SMALL FORCING STRENGTH PARAMETER 

This section considers the limiting regime where the forcing strength < 1. This 
limit is expected to be applicable to the reheating problem after an inflationary epoch 
(Appendix B). The reheating phase takes place over many oscillations of the inflaton field 
and hence many cycles of the corresponding Hill's equation. As a result, the fluctuations 
(given by the magnitude of the qk) must be relatively small. 

In general, when the forcing parameter q^ is small, solutions to Hill's equation tend 
to be stable in the classical regime, i.e., where the parameters do not vary from cycle to 
cycle. However, variations in the parameters (Xk, qk) allow for unstable solutions, even 
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Figure 1: Degree of validity for the approximations to the growth rates for Hill's equation in 
the delta function limit with large forcing parameters. The lower curve shows the difference 
between full growth rate calculated from matrix multiplication and that calculated from 
Theorem 2.1 (equation [9]). The upper curve shows the difference between the full growth 
rate and that calculated from the more extreme approximation of Corollary 2.1 (equation 
[20]). Here, the values of A are fixed and the values of qk fluctuate according to qk = 
(1 + £,k)qo, where is a uniform random variable and the constant go provides a measure 
of the fluctuation amplitude (where (qk) = 3go/2). 
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if the growth rate would vanish in the absence of fluctuations. For the case of delta 
function barriers, the parameter qk/y/Xk must be small for stability, so that classically 
stable solutions can also arise in the limit of large oscillation frequency X k (see equation 
[5]). Here we find the growth rate for a random Hill's equation in the limit of small forcing 
strength for the case of delta function barriers: 

Theorem 3.1: Consider a random Hill's equation (TjQ) with a delta function barrier so 
that Q{t) is given by equation (T5]). In the limit of small q k <C 1, fixed A, and symmetric 
variations in the q k , the growth rate approaches the form: 



where the angular brackets denote expectation values. This form is valid for all A > 
except for narrow bands of stability centered on square integer values A = n 2 (n G Z). 
The growth rate vanishes at these values of A and the width S A of the bands is given by 



Proof: In the limit of delta function barriers, the principal solutions are given by equation 
([5]). For the case of small forcing parameters < 1, the h k are less than unity except 
for the narrow zones of parameter space defined by the condition A = n 2 and by equation 
( 123]) . As a result, we can rewrite the matrix elements hk in the form 



7 = log 1 + (g fc 2 )/8A 



(22) 



5\ « 2q k /n. 



(23) 



h k = cos 9 k • 



(24) 



The transformation matrix of equation ([3]) can be written in the form 



M k 



cos Ok 
.sin 9 k /L k 



-L k sin 9 k 

cos Ok 



(25) 



where the parameter Lk is defined by 



sin 6 k sin 1 + (q k / VX) cot <p - {q\j 4A) sin 2 <p 



(26) 



9k VAsin<^ + (q k /2)(l + coscp) 

To leading order in q k -C 1, L k — Lq — —1/y/X, where L is a constant 
L k = Lq(1 + T] k ), the perturbations can be written in the form 



If we write 



sin 6k 



- 1 




+ 0{q 2 ) , 



(27) 



sin if + (g fc /2v / A)(l + cos ip) 



where the second equality defines the leading order expression. 
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As shown below, the product (^sin^) appears in the expression for the growth rate 
and is thus the quantity of interest. Since rj k is first order in q k , we can use the leading order 
expression for sin 6^ to evaluate the product. For small gt < 1 trigonometric identities 
imply the following transformation between the angle ip k and the angle Q k . 

As a result, sin 6^ = simp + 0(q k ), so that (rj k sm9 k ) = —q k /(2\/ r X) to leading order. 
Next we expand the transformation matrix of equation (125]) into two parts, 

M k = M k(O k ;L )+M lk , (29) 

where the first term is an elliptical rotation matrix with constant length parameter L and 
where the second term has the form 



Mi k = -r] k sm6 k 



L 
1/[L (1 + Vk )] 



(30) 



Note that the first term in equation ff29l is stable under matrix multiplication. Notice also 
that the second term A4\ k , as written, includes the full correction (with no approximations). 

As shown in the following analysis, the first non-vanishing contribution to the growth 
rate is second order in 7] k . As a result, we expand the product of N matrices A4 k , 

M k N) = (M ok + M lk f , (31) 

including all terms to second order in the matrix Ai\ k , 

N 

M W) = ^S ) + E^+EG£- ( 32 ) 

fc=i k,e 

The first sum includes partial product matrices of the form 



fc-i 




?k =\ II M j0 \ M lk m M io \ = £o(a k ; L Q ) M lk £ (b k ; L Q ) . (33) 
[j=k+i J [j=i J 

In the second equality we have evaluated the products using the properties of the elliptical 
rotation matrices, denoted here as Sq, and we have defined the composite angles 

N (fc-1) 
j=k+l j=l 

The second sum in the expansion of equation (|32|) involves partial product matrices with 
the form 

Qki={ II M, tJ } Mu- { II M (>j } M n { IT M (>j } . m) 
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This second sum includes all possible products of the above form, i.e., all possible locations 
of the two matrices that are of type M.\k rather than Aiok- By construction, each matrix 
Qke contains two factors of the random variable so that Qke oc rjkrje, where the rjk and r/e 
are independent realizations and hence are uncorrelated. The matrix elements from the 
third term in equation (132]) . the sum over the Q^ e , must thus vanish in the limit N — > oo 
(see also [3]). As a result, we only need to consider the contribution from the first sum in 
equation ( 1321) . In this sum, the first order terms, those proportional to rjk, will also vanish 
in the limit N — > oo. We thus need to include the second order terms in the first sum. 
Using the result of equation fl33l and expanding to second order in r)j, we thus obtain 



N N 



fc=i 



k=l 



— sin ak cos bk L sin ak sin bk 
(l/I/o) cosafc cos bk — cosafcSin& fc 



(36) 



After using this result in the expansion of equation ( l32|) . and evaluating the product -M,^, 
the eigenvalue A for the full product matrix after N steps is given by 



N 



A 2 — 2A cos On + 1 + E {jll s i n ®k [(A — cos 6^) sin ak + sin On cos ak] | = , (37) 

fc=i 

where we have defined ak = ak + bk (see equation [53]). and where 

(38) 

The zeroth order contribution to the eigenvalue is given by 



N 



An = cos On ± i sin 



>N 



(39) 



and the leading order correction is given by 



±i N 

A 2 = — E {vt sin °k (cos a k ± i sin a k ) } . 



(40) 



fc=i 



The magnitude of the full eigenvalue, A = Ao + A 2 , is then given by 

1 N 1 
|A| = 1 + r E Vl sin2 k =l + -N( V 2 k sin 2 Ok) , 
fc=i 



(41) 



where the second equality is valid in the limit A^ — > oo. To leading order in rjk, this 
expression can be rewritten in the form 



|A| 



1 

I + -( V 2 k sin 2 e k ) 



(42) 
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The corresponding growth rate thus becomes 

(43) 



N— >00 iV 



l + ^sin 2 ^) 



Using equation (1271) to determine r] k , we obtain the expression claimed in equation (122]) . 

To prove the second part of this theorem, we note that the expansion of equation ff27j) 
is no longer valid when the second term in the denominator of equation ( 1261) dominates 
the first. The condition for the expansion to fail can then be written in the form 



q k cos 2 0/2) 



\/A siny? 



(44) 



The left hand side of this equation blows up when sunp = 0, which occurs when ip = rnr 
and n is an integer; equivalently, this singularity occurs when y/X is an integer (and A is a 
square integer). Near these square integer values of A of interest, we can write 

(f = V%t = mr + 5ip , (45) 

where the second equality defines 5(p. Combining the above two results implies that 

-ISt/jL [i + (_i)« cos(^)] ~ 6<p . (46) 
nir + dip 

For n even, we thus obtain 5<p ~ gfc/n to leading order. By definition, 5X = 2n(5ip)/7t, so 
the width of the interval where equation (1221) fails is given by 5X ~ 2q k /Tt, in agreement 
with equation ( 1231) . This derivation applies to even integers n. For the case of odd n, one 
can derive the analogous result. • 

To illustrate this set of results, we present the following numerical experiment: The 
natural oscillation frequency, as set by the parameter A&, is fixed at a constant value. The 
forcing strength is then allowed to vary according to the ansatz 

Qk = qo€k , (47) 

where go is a fixed amplitude and is a uniformly distributed random variable with 
— 1 < ^ < 1. In the absence of the fluctuations, Hill's equation would have bands of 
stability and bands of instability in the (X-q) plane of parameters [1,2,19]. However, with 
the cycle to cycle variations of the forcing strength given by equation (147j) . the bands of 
stability essentially disappear. Figure [2] shows the growth rate plotted as a function of A 
for a collection of amplitudes go (where the amplitudes are equally spaced logarithmically, 
so that go = 10/2 for i = 4,5,6,7,8). Even for extremely small values of go (and hence 
correspondingly small qk), the growth rates are nonzero. The growth rate does vanish for 
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Figure 2: Growth rates for Hill's equation in the delta function limit for fixed values of A 
and fluctuating values of q^. The five curves shown here correspond to five values of the 
fluctuation amplitude go, where the qu = go 6^ where is a uniformly distributed random 
variable — 1 < < 1. For the five curves shown, the amplitudes are given by go — 10/2 f 
for i = 4, 5, 6, 7, 8. The dashed lines show the limiting form for the growth rate from 
equation (1221 . 
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Figure 3: Comparison of growth rates for Hill's equation in the delta function limit for 
fixed A and fluctuating values of the forcing strength q^. The solid curve corresponds to 
the exact result from matrix multiplication using a fluctuation amplitude go — 2.5. The 
dot-dashed curve shows the approximation developed in this section (equation [22]). The 
dashed curve shows the growth rate that results from an average of the growth rates 
for individual cycles (the asymptotic growth rate). Finally, the dotted curve shows the 
growth rate resulting from a constant value of the forcing strength q = go/2. 
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particular values of the frequency parameter A, where A = n 2 and n is an integer. At these 
particular frequencies, sin# = and hk = cos# = ±1 = constant. Notice also that in 
the limit go 1, this incarnation of the random Hill's equation is equivalent to a simple 
harmonic oscillator with frequency A and perturbative noise; this result - that noise leads 
to instability or can speed up instability - has analogs in previous work [9]. 

For the particular choice of qk used in Figure EJ (q%) = q 2 /3. The dashed curves in 
Figure [2] show the approximation of equation ( 1221) for the growth rate. In this case, the 
approximate form actually gives more accurate results than direct matrix multiplication 
(solid curves) due to incomplete sampling in the latter. However, the expression of equation 
(I22p does not account for the vanishing of the growth rate for particular values of A. 

As another illustration of how random variables change the landscape of parameter 
space, Figure [3] shows the growth rates as a function of (fixed) A for several cases. The 
solid curve shows the result from full matrix multiplication, with the same sampling of 
qk as used in Figure El with value go = 2.5. The dot-dashed curve shows the prediction 
from the approximation of equation ( |22l) . The approximation works well except near the 
integer square values of the frequency A where the growth rate vanishes. However, since 
the qk are of order unity, rather than fully in the regime ^Cl, the assumptions for the 
validity of equation ( 1221) are not completely satisfied. As a result, the small amplitude 
oscillations of the approximate result (dot-dashed curve) about the true growth rate (solid 
curve) are real. In Ref. [2] we defined the asymptotic growth rate 7^ to be the growth 
rate for a random Hill's equation resulting from an appropriate average of the growth rates 
for the individual cycles; this quantity is plotted as the dashed curve in Figure [3] For this 
set of qk values, the amplitude is g , so the mean of the forcing parameter magnitude is 
go/2; for comparison, the dotted curve shows the growth rates for the classical problem 
(no random variables) with q = go/2. The dotted curves thus delineate the regions of 
stability and instability that characterize the parameter space of Hill's equation. Note 
that the asymptotic growth rate 7^ is sometimes larger and sometimes smaller than the 
true growth rate. As a general rule, one finds 700 > 7 in or near the portions of parameter 
space for which the classical problem (fixed g) is unstable; for regimes in which the systems 
is classically stable, however, the opposite holds so that 7 > 7oo. 

IV. FOKKER-PLANCK APPROACH 

For the case of constant oscillation frequency parameter A and sufficiently small forcing 
strengths qk, Hill's equation in the limit of delta function barriers can be described by a 
Fokker-Planck equation of conventional form. Following standard methods, we define the 
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velocity V and diffusion constant D according to 

V = % and D = (48) 

where 7r is the period of the forcing intervals. The Fokker-Planck equation for the evolution 
of the distribution P(y, V, t) of phase space variables thus becomes 

dP „dP . dP D 2 d 2 P 

ln +v ^- Xy dv = J y dv*- (49) 

Notice that including variations in the or working in the regime of large forcing pa- 
rameters qk (e.g., the highly unstable limit), would require additional terms in equation 

(USD. 

In order to reduce the complexity of equation (jUJ), it would be useful to average over 
one of the independent variables. In this case, however, both V and y are on a nearly equal 
footing. In this problem, the system acts like a simple harmonic oscillator, except at the 
delta function barriers where its energy jumps due to the forcing. We thus transform into a 
type of "polar coordinates" [20] in which the energy plays the role of the radial coordinate; 
specifically we define 

E = ^{v 2 + Xy 2 ) and V = tan" 1 (y/Xy/v) , (50) 
where the corresponding inverse transformation takes the form 



V = 2VE cos ip and V\y = 2\ r Esinip . (51) 

In terms of the new variables (E,ijj), the Fokker-Planck equation becomes 
OP r-dP 2D 2 , f dP ATn 2 , dP 2 

~dt + 9^" = ~ 8111 \dE + cos +a& 

sm2ipdP sm 2 i/jdP 2 , d 2 P ) 

+ + —w - 2sm2 ^^ j • (52) 

If we now average over the angular variable ip, the partial derivative terms with respect to 
ip vanish, and the Fokker-Planck equation simplifies to the form 

dP 2D f. 2 ,,dP 2 , 2 \ dp 2 1 D f BP „ 2 dP 2 ) , , 

Next we change variables again, by defining 

fx = \ogE. (54) 
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Note that we must take the logarithm of a dimensionless quantity. Given the form of 
equation (|53|) . however, it is straightforward to introduce a dimensionless energy E = 
E/Eq before changing to the logarithmic form. The resulting Fokker-Planck equation thus 
becomes an ordinary diffusion equation 



OP D dP 2 



dt A dy 2 

This diffusion equation has the normalized solution 

r,, ^ ( A 
P(fX,t) 



(55) 



AuDt 



,1/2 


" /i 2 A" 


exp 


4Dt 



(56) 



which is appropriate for the boundary condition P(fi,t = 0) = 5(/i), i.e., the system 
starts out with ji = or energy E = E . The solutions y{t) are oscillatory, but growing 
(in general). In order to extract a growth rate from this Fokker-Planck treatment of the 
problem, we first determine the expectation value of y 2 , which takes the form 

(y 2 ) = < /_ oo t)dn = - exp [Dt/ A] . (57) 

To obtain the final expression, have averaged over the angle ip, so that 2 sin 2 ?/' = 1. Next 
we assume that the amplitude \y\ of the solution can be characterized by the relationship 

|y|^(y 2 ) 1 / 2 ocexp[ 7fp t] , (58) 

where 7f p is the growth rate resulting from this Fokker-Planck approach. The resulting 
estimate for the growth rate is 

This growth rate is similar, but not identical to, that given by Theorem 3.1 for the case of 
constant frequency parameter A, delta function barriers, and in the limit of small forcing 
parameters q^. The functional dependence 7 oc q 2 /X is the same, only the numerical 
coefficient differs. In order to derive the result (1591) . however, we have averaged the Fokker- 
Planck equation ( |52l) . and this procedure can produce such a numerical difference. 

This treatment using the Fokker-Planck equation thus provides a good description of 
the problem for the case of small qk and constant A. For varying values of A&, additional 
diffusive terms must be included. For the case of large qk, however, the Fokker-Planck 
approach does not naturally reproduce the results obtained here using direct methods. 
The derivation of equation (H91) involves truncating a series of terms in powers of q^ [6], 
and such a truncation is only valid for sufficiently small forcing strengths. In the highly 
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unstable limit (large qk), one must either use a highly modified form of the Fokker-Planck 
equation or abandon this approach altogether. 



V. TRANSFORMATION FROM HILL'S EQUATION PARAMETERS TO 
MATRIX ELEMENTS 

For Hill's equations with delta function barriers, we can write the matrix elements 
hk and g k in the form given by equation (jSJ) above. With these results in hand, we can 
directly construct the relationship between the fundamental parameters (A/-, qk) appearing 
in the original differential equation ([1]) and the moments of the distributions of the matrix 
elements. 

We start by considering the case where the angle (p k is held fixed, but the forcing 
strength qk is allowed to vary. For ease of notation, we suppress the subscripts on the 
angle (p and the forcing strength q. The mean (first moment) of the matrix element hk is 
then given by 

. , . 7r . . sin <p> . , 

(h k ) = coa<p--(q)—Z. (60) 

L <p 



Similarly, the second moment takes the form 



</4> = COSV + ^<S 2 > I — J -7T — (g). (61) 



2 



/ sm (p \ cos (p sin ip 



The variance is thus given by the expression 

2 / • \ 2 

n I sin 9? \ r 2 2l 7r 



4 V if I 1 J 4 V (f 




where the second equality defines the variance of the distribution of the forcing strength 
q. As a result, this limiting case provides a simple relationship between the width of the 
distribution of forcing strength q and the width of the distribution of the resulting matrix 
element hk, i.e., 

a h = ^—a q . (63) 

Next we consider the opposite case in which the forcing strength q is held fixed and the 
angle (p varies over a range. Here we consider a range of angles so that expectation values 
are taken via the operator 

(...) = i fdip... , (64) 



o 



r 

which holds for any given quantity in the brackets. Notice that we are using a distribution 
that is uniform in the variable (p. Since (p oc \/A, this distribution is not uniform in the 
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variable A, although an analogous analysis could be done for that case. For this choice of 
distribution, the first moment of the matrix element is then given by the expression 

<fc> = =^->(r). (65) 

where Si(T) is the sine integral [1]. The second moment is given by 
2 1 | sin2r | ( nq\ 2 \ 1 a . friT ^y sin 2 r] 7rg ( 



<*» = 2 + — Hf J |r Sl < 2r > " — 1 " 2F Sl < 2r » ' < 66 > 

These (general) expressions can be simplified by choosing the angle interval to be have the 
form T = 27rm, where m is an integer. In this case, the corresponding variance reduces to 
the form 

^ = i + S^-iWr)-^y S i>(r). (67) 



2T V 2 ) y ' \2TJ 
Next we note that these results simplify further in the limit T — > oo, i.e., when we allow 
the angle tp to vary uniformly over the entire positive real line. In this limit we find 

1 \/2 
lim (h k ) = , lim (h 2 k ) = - , and lim a h = — . (68) 
i ^oo r^oo z r^oo z 

As long as the distributions of the angle ip and that of the forcing strength parameter 
q are independent, the expressions derived above for the moments can be generalized to 
include both distributions in a straightforward manner. In particular, we continue to use 
the distribution of equation (|6%|) for the angle, and the same (unspecified) distribution 
for the forcing strength as above. In this case, we must integrate over both the forcing 
strength q and the angle (p. By performing the integrals over q first, we obtain expressions 
analogous to equations (|60|) and (|6T!) for the first two moments. With these results in 
hand, we then average over the distribution of angle using equation (1641) . This procedure 
produces expressions of the forms given by equations ( 1651) and ( 1661) . with q replaced by (q) 
and with q 2 replaced by (q 2 )- For the case in which T = 2irm, the resulting expressions for 
the first two moments take the form 

(h k ) = ~&F)(q) and (h 2 k ) = \ + ^Si(2r)(g 2 ) - ^Si(2r)(g> . (69) 

The corresponding variance is thus given by 

4 = l + ^Si( 2 r) {(q 2 ) - 2 -{q)) - ^(g) 2 Si 2 (r) . (70) 

If we consider the limit where both (q) and T are large, the variance of h k has the order 

o* h = l + (a 2 /T) . (71) 
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Next we consider the analogous calculation for the matrix element g\.. For the limiting 
case in which the angle is held fixed and the forcing strength varies, the first two moments 
of the distribution are given by 



and 



i i 

(9k) = V 5 s in <P ~ o ( X + cos V) (?) 

7T 2 



11 1 

-ip 2 sin 2 (p + — (1 + cos (p) (q 2 ) (p sin <f(l + cos ip) (q) 



7T" 



7T 



The corresponding variance thus takes the form 

l(l + cos^) 2 {(g 2 )-(g) 2 } 



(72) 
(73) 

(74) 



5 4 

This result can be rewritten in a manner analogous to that found for the other principal 
solution (equation [65]). i.e., 

cr g = - (I + cos ip) a q . (75) 

For fixed forcing strength q, and a distribution of angle given by equation (JM]) . the first 
moment of the distribution becomes 



<0fc> 



1 / sin T 



7T 



r 



cos r 



sin T 



(76) 



and the second moment is given by 



(9l) 



7T 



T 2 cos(2r) 2r 2 - 1 



6 ir 

q 



8T 



sin(2r) 



+ 



3 2 sin T sin(2r) 
- H 1 — - 

2 r 4r 



+ 



7T 



„ 1 , sinT sin(2r) 

- cos r - - cos 2r + — + — i — '- 
4 v i r sr 



(77) 



These expressions are somewhat cumbersome; if we take the angular interval to be T = 
27rm, where m is an integer as before, all of the sine terms vanish and the moments simplify 
to the forms 



{9k} = - ~ - 

7T Z 



and 



1 Y 



3q 2 5q 



(9k) = -dir-7F+^r-7Z 



7T- 1 



6 4r 



4tt 



In this case, the variance is given by 



'r 2 



q 



a " 7T 2 \ 6 4r 1 J + 8 4?r ' 



(7f 



(79) 



x 6 4r 

In the limit V — > oo, the width of the distribution does not converge, but rather diverges 
linearly so that 

(80) 



hm o„ ~ — - 1 



7r6 
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As before, we can simultaneously include the distributions of angle and forcing parame- 
ter, provided that the variables are sampled in an independent manner. For the case where 
the angular interval is taken to be T = 2irm, the moments reduce to the forms 

±) + h g 2)_ JL 2 (81) 
AT) 8 W 1 4vr 7 v ; 

<?> 2 }-^>- (8 2 ) 

As a way to compare with the results found for the moments of the h\. distribution, we 
again consider the limit where both (q) and T are large. In this case, the variance of the 
gk has the order 

a 2 = 0(T 2 ) + 0(a 2 ) . (83) 

Comparing this result with equation (J7II) . we find that when the angular interval V is not 
too large, the variance of the hk and that of the are of the same order. When the interval 
T is large, however, the variance of the g^ dominates. 

Next we consider the limit of large forcing parameters qk- In this limit, the leading 
order growth rate is given by Theorem 2.1, and the corrections are determined by the 
variables Xk and <pk appearing in the matrix of equation ([7]). For the sake of definiteness, 
we define the mean q Q = (qk), with q large, and allow the qk to have large variations about 
the mean, but not so large that qk — > 0. Under these conditions, one can show that the 
variances of the variables Xk and <pk (from equation [7J) are of the order 

a 2 x = O(q 2 ) and a\ = O(q^) . (84) 

As a result, the variance is dominated by the Xk rather than by the fa. This result is 
consistent with the findings of Section II. 

VI. CONCLUSIONS AND DISCUSSION 

This paper has generalized and extended previous work concerning Hill's equations 
(CQ) that contain random forcing parameters, with a focus on the case of delta function 
barriers (equation [2]). In this formulation of the problem, both the natural oscillation 
frequency Xk and the forcing strength qk can vary from cycle to cycle. The development 
of the solutions to Hill's equation, including the growth rates for instability, are given by 
the general problem of matrix multiplication (equation [3]), where the matrix elements are 



11 1 (Y 2 

(gk) = - - -(q) and (Pfc) = — -7- - 

7T 2 7T Z y o 

The variance for this case is given by 



2 _ J_/T!_ J__ " 

Gg ~ 7T 2 I 6 if ~ 1 
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determined by the principal solutions for a given cycle. We have constructed the principal 
solutions for individual cycles using Dirac delta functions as the periodic barriers (equation 
[5]). This construction allows us to explicitly show how both the matrix elements hf. and 
gk, and the growth rates 7, depend on the distributions of the original parameters (Afc, g&) 
appearing in Hill's equation (Section V). 

In the limit of large forcing strength parameters g&, the growth rates approach the 
form 7 ~ (log |<7fc |) (see Theorem 2.1). In this limit, large qt values often lead to large 
values of the principal solutions at the end of the cycle; this result, in turn, demonstrates 
(by construction) that the highly unstable limit [2] can be realized using unremarkable 
values of the parameters. In the opposite limit of small forcing parameters, the growth 
rates approach the form 7 ~ (ql) (Theorem 3.1). Our results in this limit show how the 
fluctuations in the Hill's equation parameters act to fill in the bands of stability in the 
classic problem with fixed parameters (Figure [3]). We have also found analytic results 
for the widths of the remaining bands of stability (where the growth rates vanish). In 
this same limit (small qk), the Fokker-Planck equation provides an alternate description 
of the dynamics (Section IV) and consistent estimates for the growth rates. Finally, we 
have constructed a iterative map (Appendix C) that provides a heuristic argument for the 
general form of the growth rates in the limits of both large and small forcing strength q^. 

The results of this work can be used in a number of applications. For example, one 
motivation for considering random Hill's equations was to study orbital instabilities in 
extended mass distributions, such as dark matter halos, galactic bulges, and young embed- 
ded star clusters (see Appendix A). The results presented herein show when the orbits are 
unstable and provide estimates for the corresponding growth rates. These results, in turn, 
help explain the observed dynamical structures in these astrophysical systems. Another 
important application involves the reheating problem at the end of the inflationary epoch 
in the early universe (see Appendix B). In this context, the introduction of stochastic 
perturbations (e.g., due to quantum fluctuations) leads to the disappearance of the bands 
of stability (see Figures [2] and [3]) . As result, fluctuations enhance the effectiveness of the 
reheating process. In addition to these motivating examples, random Hill's equations arise 
in a wide variety of other physical problems [2-5,7,14,15,19-21]. 
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APPENDIX A: 

RANDOM HILL'S EQUATION FROM ASTROPHYSICAL ORBITS 

One application of Hill's equation with random forcing terms involves the study of an 
instability that affects orbits in extended mass distributions, such as dark matter halos [4]. 
In this setting, the density profile p{vo) of the halo has the general form given by 

p{w) = po ^ , (Al) 

where p is a density scale and the variable w is written in terms of the usual (x, y, z) 
coordinates through the relation 

^2 x 2 y 2 z 2 
cr b l & 

where a > b > c > 0. The density field is thus constant on ellipsoids. The function F(zu) 
is approaches unity as w — > so that the density profile approaches the form p ~ 1/zu in 
the inner limit. For this regime, one can find analytic forms for both the potential and the 
force terms [4]. 

Further, when an orbit begins in any of the three principal planes, the motion can 
be unstable to perturbations in the perpendicular direction. Consider an orbit initially 
confined to the x — z plane, with a small perturbation in the perpendicular y direction. In 
the limit \y\ <C 1, the equation of motion for the y-coordinate takes the form 

4^ + u?.y = where ul = , . (A3) 

dt 2 yy y ^c 2 x 2 + a 2 z 2 + bVx 2 + z 2 1 ; 

In this setting, the time evolution of the coordinates (x, z) is determined by the origi- 
nal orbit. Since this orbital motion is nearly periodic, the [x(t), z(t)] dependence of the 
parameter uj 2 provides a periodic forcing term. The orbit has a maximum extent (outer 
turning points) which results in a minimum value for u 2 , which in turn defines the natural 
oscillation frequency A^. The parameter oj 2 defined above can thus be written in the form 

wl = . A ' b =A fc + Q fe (t), (A4) 

y Vc 2 x 2 + a 2 z 2 + bVx 2 + z 2 K J 
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where the index k counts the number of orbit crossings, and the chaotic orbit in the 
original plane leads to different values of Xk and Qk(t) for each crossing. The shape of the 
functions Qk are nearly the same, however, so that one can write Qk(t) = qkQ(t), where 
the forcing strength parameters qk vary from cycle to cycle. These forcing strengths qk are 
determined by the inner turning points of the orbit (with appropriate weighting from the 
axis parameters [a,b, c]). Given the expansion of equation (1A4I) . the equation of motion 
( 1A3I) for the perpendicular coordinate becomes a random Hill's equation, with the form of 
equation ([!]), as studied herein. 



APPENDIX B: 

RANDOM HILL'S EQUATION FROM REHEATING IN INFLATION 

In the inflationary universe paradigm [12], the accelerated expansion of the universe is 
(usually) driven by the vacuum energy associated with a scalar field <p (often called the 
inflaton). During the phase of accelerated expansion, the energy density of the universe 
itself decreases exponentially and the cosmos becomes relentlessly empty. This epoch is 
thought to take place when the universe is extremely young, with typical time scales of 
~ 1(T 36 sec. 

In order for the inflationary epoch to solve the cosmological issues it was designed to 
alleviate, the end of inflation must involve a mechanism to fill the universe with energy 
(e.g., see the review in Ref. [17]). This process is called reheating. During the epoch of 
reheating, the equation of motion for the inflaton field displays oscillatory behavior about 
the minimum of its potential. Further, in order for the universe to become filled with 
energy (reheat), the inflaton field <p must couple to matter or radiation fields. One simple 
type of interaction that is often considered uses an coupling term in the Lagrangian of the 
form 

Ant = 5 , ^X 2 , (Bl) 

where x is a second scalar field that represents matter (radiation) and where the coupling 
constant g sets the strength of the interaction. The field \ is generally expanded in terms 
of its Fourier modes Xk since these quantities evolve independently. The resulting equation 
of motion for the matter field modes Xk takes the form 

d 2 Xk 



dt 2 



+ 



oo 2 k +p(t) + q(t) X k = 0, (B2) 



where p(t) is a periodic function (given by the oscillatory behavior of the inflaton field) 
and q(t) is a noise term that provides perturbations to the driving term p(t) [21]. Note 
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that the index k refers here to the Fourier mode, although the forcing terms do vary from 
cycle to cycle. In the absence of fluctuations, the matter field modes Xk thus obey a type of 
Hill's equation, which is subject to parametric instability [15,16]. The noise perturbations 
convert the equation into a random Hill's equation [15,16,21], of the type studied herein. 
This type of equation was solved numerically using WKB methods [16], thereby finding 
the relevant physical solutions; nonetheless, the formulation of this paper can be applied 
to this class of reheating problems, and more general results can be obtained. 



APPENDIX C: AN ITERATIVE MAP 



As shown in the text, the growth rates for Hill's equation depend on the forcing strength 
q k according to 7 ~ (qf) in the limit of small symmetric q k , and 7 ~ (log \qk\) in the limit 
of large q k . These results hold both for the particular case of delta function barriers 
(considered here), and for the general problem [3]. In this Appendix, we construct a 
heuristic argument that reproduces these forms for the growth rate in the two limits. This 
treatment is highly approximate, by design, but allows for a simple interpretation of our 
previously obtained results. 

Given the form of Hill's equation in the delta function limit, the jump condition across 
the barrier takes the form 

dy 



dy 
dt 



dt 



-QkV, 



(CI) 



where all of the functions are evaluated at the barrier. If we define V = dy/dt, and relabel 
the functions with an index k + 1 on the far side of the barrier, and an index k on the near 
side, we obtain an iterative map of the form 



k+l 



v k 



Vk 

qt v t 



N 



Voll 

k=l 



Vk 

qt v k 



(C2) 



where we have continued the iteration back to the initial step to obtain the second equality. 
The growth rate 7 for this map can then be defined according to 



1 



N 



7 = lim — lo 



k=l 



1 - qk 



Vk_ 

Vk 



(C3) 



Given the form of Hill's equation away from the delta function barrier, the solutions are 
oscillatory with frequency y/Xk, so that the function y k and the velocity are related via 

Vk 1 
Vk 



FUXkt), 



(C4) 
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where the function F depends on the angle \/Xk t. If we use the ansatz implied by equation 
( 1C4|) . the growth rate takes the form 



7 



lim -J- 

N-*oo N 



N 



k=l 



1 - 



Qk 



(C5) 



In the limit of large forcing strength ^> 1, the growth rate reduces to the form 



iV 



7 



lim — 



fe=i 



log 



Qk 



(log |/lfc|) 



(C6) 



where we have ignored the function F of the angle in obtaining the final approximate forms. 
This argument reduces the problem to a single (approximate) iterated jump condition, but 
still reproduces the proper dependence of the growth rate for the highly unstable limit 
(7 ~ (log \ hk\))- We note that ignoring the angular function is not valid when F — ► 0. As 
shown above, this problem allows for narrow bands of stability where the growth rate can 
vanish even when the forcing strength is large. The presence of stable behavior (7 — > 0) 
can thus be accounted for through this heuristic argument (by allowing F — > 0). 

In the opposite limit of small forcing strength <C 1, we can expand the logarithmic 
function in the expression for the growth rate to obtain 



N 



7 



lim — 

,r .„ AT ^-^ 



k=l 



Qk 



■ F + 



qI 

2X k 



(C7) 



For symmetric fluctuations, the first term vanishes in the limit, so that the quadratic term 
provides the leading order contribution to the growth rate. Ignoring the angular function 
F as above, the growth rate becomes 



7 



Qk/h 



(C£ 



In this case, the iterated jump condition argument reproduces the proper dependence of 
the growth rate for the limit of symmetric and weak forcing (compare with Theorem 2.1). 
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